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E. MUKHIN* AND V. TARASOV* 

Abstract. We study polynomials p{x) satisfying a differential equation of the 
form p" — h'p' + Hp — 0, where h = x 3 /3 + ax. We prove a conjecture of A. 
Eremenko and A. Gabrielov. 



1. Introduction 
In this note we consider the equations of the form 
(1.1) p"{x) - (x 2 + a)p + H(x)p = 

which have a polynomial solution. Such equations appear in the study of the ele- 
mentary eigenfunctions of the Schrodinger equation with quartic potential, see |EG| . 

We consider the corresponding local system. The cohomology of the system is two- 
dimensional. Our main result is the proof of a conjecture of [EGj . which describes 
the cohomology class of the polynomial p 2 (— x). 

The problem of computing the cohomology classes is formulated in an algebraic 
setting, see Section [31 but we use complex-analytic tools to solve it. Our main 
insight comes from the consideration of the bispectral dual equation to (11.11) . see 
equation (15.11) . 

To prove the wanted equality of two constants we interpret them as values at 
zero of two a priori different solutions of equation (15.11) and then show that the two 
solutions actually are the same comparing their asymptotics via steepest descent 
method. 

The paper is written as follows. In Section [2] we discuss the local system associated 
to (II. ip . We proceed to describe an explicit basis in the cohomology in Section [31 
In Section H] we exhibit polynomials which are homologoues to a constant multiple 
of the first basis element proving in particular Conjecture 1 from [EGj . We discuss 
the bispectral dual equation in Section Section [6] is devoted to the elementary 
computation with the characteristic equation of the linear operator corresponding to 
(11.11) . Finally, we prove our main results Theorem 17.11 and Corollary 17.21 in Section 

m 
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Figure 1. Contours jj and planes Hj. 



2. Elementary remarks 



Fix a G C and let 



h(x) 



x 



— + ax G C[x] 



Denote by the prime the operator of differentiation with respect to variable x and 
define a linear map on rational functions of x: 



D : C{x) ->• C(x) 



i — y q'(x) + h'(x)q(x). 



d_ 

dx 



The map .D is inherited from the derivative map 
Let C C C(x) be the image of D. Let 

R = {q{x) G C(x) | Resg(a;)e^ ) 



[x e 



0} 



ft(a;) 



->■ C(x)e 



ft (a;) 



be the subspace of rational function which have no residues after multiplication by 
the exponential of h(x). We have C C R. 

For qi(x), q2(x) G C(x), we write qi(x) ~ 92^) if and only if qi(x) — q%{x) G C . 

Let Jj(t), j = 0, 1, 2, be any smooth curves in complex plane such that 



lim arg( 7j (t)) = 7r(l/3 + 2j/3), 

t—¥— 00 



lim arg( 7i (t)) = vr(l + 2j/3) 



t— >oo 



and lim t _s. ±00 |7j(t)| = c>o, see Figure [TJ 

Define functionals lj G R*, j — 0, 1, 2, by the formula 

^(x)) = / g(x)e h(a; W. 



7 3 - 



Here we chose the contour jj so that it does not go through possible poles of q(x). 
Clearly, the functionals ^ are well-defined and we have l\ + I2 + I3 — 0. 
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Proposition 2.1. We have dim R/C = 2. Moreover, for q(x) G R we have q(x) ~ 
if and only if lj(q(x)) = 0, j = 1, 2. 

Proof. Let q(x) G R. Write q(x) as a sum of simple fractions. If we have a term 
l/(x — z) k , for some z G C, then we subtract D(l/(x — z) k ~ l ) and decrease the order 
of the pole modulo C. Note, that since q(x) G R, k = 1 is impossible. Therefore 
there exists a polynomial f(x), such that q(x) ~ /(^)- Now if deg/(x) = n > 1 
we subtract D(x n ~ 2 ) and reduce the degree of f(x) modulo R. It follows that 
q(x) ~ ax + b for some choice a, b G C. 

Since |e fe ( 7; >'^| — > as t — >■ ±oo, we have /j(C) = 0. Therefore to finish the 
proof of the proposition, it is sufficient to show that G R* are linearly inde- 
pendent functionals. Thus, it is sufficient to show that det(/^ x k ~ x e h ^ dz\^ k=x ^ = 
0. But this determinant is non-zero because it equals the Wronski determinant 
W(4>i(u), 4>2(u)) where 4>j{u) = e h ^ +ux dx are fundamental solutions of the Airy 
equation f"(u) + (u + a)f(u) = 0. □ 

We remark that one can replace the cubic odd polynomial h(x) with an arbitrary 
polynomial of degree k and similarly define k functionals and prove a generalization 
of Proposition 12.11 with dim R/C = k — 1. 

Let p(x) G C[x] be a polynomial of degree n with simple roots only. Let 

R(p(x)) = R n ^r~y c( P (x)) = chd . 

p z [x) \P{x) J 

Then clearly C(p(x)) C R(p(x)). If qi(x),q 2 (x) G R(p(x)) C R then clearly q\(x) ~ 
?2 (a^) if and only if gi(x) — 92(2^) G C(p(a;)). 

Lemma 2.2. VFe have dim R(p(x))/C(p(x)) = 2. Moreover, for q(x) G R(p(x)) we 
have q(x) G R(p(x)) if and only if k(q(x)) = 0, i = 1,2. 

Proof. The proof is similar to the proof of Proposition 12.11 □ 



3. When p(x) is a wave polynomial 
We start with the following lemma. 

Lemma 3.1. Let p(x) G C[x] be a polynomial with simple roots only. We have 
l/p 2 (x) G R if and only if there exists b G C such that p(x) is a solution of the 
equation 

(3.1) y"(x) - h'(x)y'{x) + {nx + b)y(x) = 0. 

Proof. Let z G C be such that p(z) = 0. We compute 

Res _ ^1 = ]im / e^)(x-^) 2 V = ch(z) h'(z)p'(z)-p"(z) 



P 2 (x) x^-z \ p 2 (z) J (P'(z)) 3 

by applying the L'Hopital rule. The lemma follows. □ 
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We call a polynomial p(x) satisfying ( 13. ip the wave polynomial of degree n. It is 
known that for each n G Z>i, a G C, there exists at least one wave polynomial of 
degree n, see also Section [6] below. 

Note that all roots of all non-zero wave polynomials are simple. 

Lemma 3.2. Let p(x) be a wave polynomial. Let f(x) be the polynomial such that 
f'(x) = p(x), /(0) = 0. Then f(x)/p 2 (x) G R(p(x)) and for any q(x) G R(p(x)), 
there exist unique c, d G C such that q{x) ~ (c + df(x))/p 2 (x). 

Proof. The claim that f(x)/p 2 (x) G -R is checked similarly to the proof of Lemma 
13.11 It follows that for any c, d G C, we have (c + df (x)) / \p(x)) 2 G R. 

If degp(x) = n then deg/(x) = n + 1. If g(x)/p 2 (x) G C then degg(x) > n + 2. 
The lemma follows from Lemma 12.21 □ 

In what follows we study the constants c, d for q(x) of the form p(— x)r(x), where 
r(x) is a polynomial of degree at most n. In particular, we prove Conjecture 1 and 
formula (18) from |EG] describing the constants c, d for q(x) = p 2 (—x). 



Let p(x) be a wave polynomial of degree n. 

Theorem 4.1. For any polynomial r(x) with degr(x) < n, there exists c G C such 
that r(x)p(—x) ~ c/p 2 (x). 

Proof. For j = 0,1,2, consider 



where the integration is taken over the contour / -fj^ x (t), t G (— oo,0], such that 
7j>(0) = x and lim^-oo arg7 j>x (t) = 7r(l/3 + 2/3j), lim^-co \jj, x (t)\ = oo. Clearly 
Hj(x) are holomorphic functions satisfying (13. ip . Let Hj C C be the half-planes 
given by 



4. When the constant d is zero 



(4.1) 




Hj = {z G C | tt(— 1/6 + 2j/3) < arg^ < tt(5/6 + 2j'/3)}, 
see Figure [TJ We have the following asymptotics: 

e h(x) 

It implies the following connection formulas 

Vj+i( x ) = VA X ) ~ J jP( x ), J = °> 1> 2 > 



where 




and y 3 (:r) 



2/0(2;), J 3 = J - 
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By Proposition [2J] and Lemma [2.21 it is sufficient to prove 

f x k p(—x)e h ^ dx J a x k p(— x)e h ^ dx ~ ^' fc — 0, 1, . . . , n. 

Make the change of variables x — >■ — x in the integrals and using the connection 
formulas we obtain that the determinant up to a sign is equal to 

J 2 / x k p(x)e~ h{x) dx-Ji [ x k p(x)e- h{x) dx 

x k (y - y 2 )e- h{x) dx- f x k (y 2 - yi )e~ h ^ dx 
71 J 72 

= / x k y e~ h{x) dx + I x k yi e~ h{x) dx + [ x k y 2 e- h{x) dx. 

J*n ^72 Jjo 

Here jj are the contours given by 7j(t) = — Jj{t) for all tel. Note that 7 J+ i C 
j = 0, 1, 2, and therefore the contour of integration in each of the last three integrals 
can be sent to infinity inside of Hj. It follows that each of the three integrals is zero 
due to the asymptotics of Vj(x). □ 

5. BlSPECTRAL DUAL EQUATION 

Motivated by |MTV] we consider the bispectral dual equation to (13. ip : 

(5.1) ug(u) — ng(u) — (u 2 — au + b)g(u) = 0, 

where the dot denotes the derivative with respect to the variable u. 

Equation (15. ip is obtained from (13. ip by formal replacing the operator of mul- 
tiplication by x with the operator d/du and the operator d/dx with operator of 
multiplication by u and placing all derivatives d/du to the right of the operators of 
multiplication by u. 

The solution sets of bispectral dual operators are often related by suitable trans- 
forms. We describe such transforms for the case of bispectral dual operators (13. ip 
and fl5TTj) . 

Lemma 5.1. Let p(x) be a polynomial solution of (13. ip . Then for j = 0,1,2, the 
integral 



gf{u) 



p(-x)e h{x) - ux dx 



is well-defined and gf\u) is a holomorphic solution of (15. ip . 

Proof. The integral is well-defined since e h( - x ^ is decaying along 7^. We twice use the 
integration by parts to compute 



ugf^ — ngf^ — (u 2 — au + b)gf^ = u{e h ux )' p{— x) — (nx — b)p(—x)dx 

'a 



(p(-x))'e h (e ux )' - {nx - b)p(-x)dx = / (p(-x)e h )' - (nx - b)p(-x)dx = 0. 



7j J Jj 
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□ 

Lemma 5.2. Let yj(x) be a solution of (13. ip given by (14. ip . Then the integral 

g f{u) =u n+l [ y,j(x)e- ux dx 



is convergent for u E C such that Re(i/e l7r ( 1+2j ' 3 )) < and g(u) is a solution of 

flHUD- 

Proof. The integral converges as t — > — oo on 7j(i) since is decaying and for 

t — t- oo since e~ ux is decaying. 

Similarly to the proof of Lemma 15.11 we use twice the integration by parts and 
obtain 

-^-riugf - ngf - (u 2 - au + \)gf) = [ (uh! - (n + 2)x - u 2 - X) yj e~ ux dx 

= I -(-^ + h' Vj ) (e- ux )' - ((n + 2)x + X) yj dx = 0. 

□ 

Using integration by parts the function gf\u) can be rewritten as follows: 
(5.2) 

where p( r )(;r) denotes the r-th derivative of In particular, the integral on the 

right hand side of (15 .2D converges for all values of u £ C and the function #j (it) is 
holomorphic in C. 

Proposition 5.3. For j = 0, 1,2, we have gf\u) = (—l) n gf\u). 

Proof. We compute the asymptotics of gj using the steepest descend method sim- 
ilarly to the computation of asymptotics of the Airy functions, see [5]. We obtain 
forj = 0,l,2, 

g f(u) = t(-l) n+ ^u n ^ 4 e-l u3/2+aul/ \l + 0(1)) 

as \u\ — > oo, with argit fixed such that 

7r/3 <argu < 7n/3 (j = 0), 

— 7n/3 < argix < — 7r/3 (j = 1), 

— 7r <argix < 7r (j = 2). 

Here j' = for j = 2 and j' = 1 for j = 0, 1. 

Similarly, using (I5.2p . we conclude that the function gf\u) has the asymptotics 

different from that of gf\u) only by a factor of (—1)™ and since there is a unique 
solution of (15. ip with such asymptotics, the proposition follows. □ 
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Lemma 5.4. Let g(u) be a solution of (15.11) holomorphic at u = 0. Then 

g(iije ux 
p{x) = Res u=0 — 

is a polynomial solution of ( 13.11) . 

Proof. We again use twice the integration by parts: 

... ,, . . , . , f (u 2 — au — x 2 u + nx + b)qe ux , 
2ni (p" - tip' + (nx + b)p) = / - — x - du 

J\u\=e U 

((u 2 — au + b)g + ng — ug)e ux 

~~ LI LL — u. 

7 ,n+l 

|it|=e a 



6. Some linear algebra 



□ 



Let V = C n+1 be the vector space with a scalar product. We fix an orthonormal 
basis {e , . . . , e n } in V. For v e V we denote u» = v i ■ e« the coordinates of v. For 
an operator A 6 End(K) we denote Ay = • Aej the matrix coefficients of A. We 
also denote A the adjoint operator of A. We have AA = AA = (det A) I. 

Let A : V — > V be a linear operator with eigenvalue —b. Let v and t> * be the 
corresponding eigenvectors of A and A*. We have (A + &)t> = 0, and (A* + 6)t>* = 0. 

Lemma 6.1. Assume that VjV* k ^ 0. Then 



d 



v ■ v 



det(A + A) \ x=b = ——(A + b) jk . 



dX J VjV* k 

Proof. Since —b is an eigenvalue, there exists a G C such that (A + b) s i = av s Vi for 
all s, I = 0,1, ... , n. Thus 

d \, . -r-r A . (^+6)^ , 

— v -v . 



— det (A + A) J | A =6 = tr aT& = « ^ 



,=0 V ^ 

□ 

We apply Lemma I6A1 to the case V = C n [x] the space of polynomials of degree 
at most n and 

(6.1) A = (d/dx) 2 - ti(x)(d/dx) + nx. 

Clearly A is a linear operator which preserves V. We choose the basis of V as 
follows: let 

efc(x) = x k /k\, k = 0,...,n. 
We set e_i(x) = e_ 2 (x) = 0. Then 

(6.2) Ae k = e k -2 - ae k -i + (n - k)(k + l)e k+1 , k = 0,...,n. 

Clearly, there exists a wave polynomial p(x) of operator (13.11) if and only if —b is an 
eigenvalue of A. Moreover, in such a case the rank of D + b is n, p(x) is unique up 
to a multiplicative constant and the degree of p(x) is exactly n. 
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Let p(x) = J2™=oPs e s(x), Ps G C, be a wave polynomial: (A + c)p(x) = 0. Then, 
clearly, p*(x) = YT s =oPn-sG s {x) satisfies (A* + c)p*(x) = 0. 
Using Lemma [6.11 with j = n, k = 0, we have 

(6.3) QL det (A + A)) | A=C = (-i)«M! |>Pn- s . 

7. The constant c 

We are now ready to compute the constants c. 

Theorem 7.1. Let p(x) = Yl™=oPs e s{x) be a monic wave polynomial of degree n. 
Then for k — 0, . . . , n we have 

(-l) n Pn-k 



e k (-x)p(-x 

EI21 th 

and for j = 0, 1, 2 



p 2 (x) 

Proof. By Lemma 12 .2\ there exists a constant Ck such that e k {— x)p{— x) ~ Ck/p 2 (x) 



JjC k = / e k (-x)p(-x)dx. 

Choose any j G {0, 1,2} and set gj{u) = gf\u) = (— l) n gf\u), see Proposition 15.31 
Using the presentation gj(u) = g^(u), we obtain 

J iCfc = y e fe (-x)p(-x)e dx = — , 

where g*f\u) denotes the fc-th derivative of gj(u). 

Expanding gj(u) in the Taylor series at u = and using Lemma [5.41 we compute 

00 c/ s) (0)M s ~ n ~ 1 00 g {s \o) 

ap{x) = Res u=0 — j e ux = ^ -?— } — e n _ s (x), 

s=0 S ' s=0 S ' 

Since p(x) is monic, we have p n = n\ and the constant a is given by gj(0)/n\. It 
follows that 

r _ sf(o) _ ^(0) 

fc tit it Pn—k- 

k'.Jj nlJj 

Finally, using that gj(u) = {—l) n gj\u) and equation (I5.2p . we obtain 

9j (0) = (-lTnUj. 

The theorem follows. □ 
Corollary 7.2. Let p(x) be a monic wave polynomial of degree n. Then 



\\=b, 



where A is given by (16.11) or (16. 2p . 
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Proof. By Theorem I7.1[ we have 

Ew.H]pH)~(-i)"^f 



p 2 { — x) 



,s=0 

The corollary now follows from (16. 3p . □ 

Corollary 7.3. Conjecture 1 and formula (18) in |EG] is true. □ 

Proof. Conjecture 1 and formula (18) in |EG] follow from Theorem 14 . 1 1 and Corollary 
17.21 respectively after the change of variables: 

x = (3z E g, p{x) = Peg(z)(3 u , (3a = 2b EG , f3 2 b = 2a EG , 

where (3 3 = —2 and we denote the objects from |EG] by the same letters as there 
but with the index EG to distinguish from the notation used in this note. □ 
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